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A ship in steering is considered as a physical mechanism that is 
۱۳ ىن‎ by a rudder movement to produce a response. “Mimphasig ds 
laid upon the relation between the forcing and the response, leaving 
aside any detailed consideration of the forces concerned. 

Course-Keeping with Automatic Control techniques is studied 
following this concept. Computer programs are developed to simulate 
different conditions. Interpretation of the results is made to evaluate 


the different methods used. 
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I. INTRODUCTION 


Problems that justify the use of Automatic Control for ships: 

The maneuverability of a ship is determined by the physical 
properties of the ship in its environment and the physiological and 
psychological properties of its navigator. 

A successful mission of a vehicle in its operating environment is 
L function of its handling as wellas its design. The relationship 
between helmsmanship and the steering system has been discussed in 
detail by Abkowitz [1], Brard [2], Davidson [3], Nomoto [4], among 
Simers, and the conclusions indicate that the hydrodynamic design 
factors alone can not significantly improve the handling qualities of 
ships. 

The past decade will be known as the decade of the Super- Tankers. 
Ships were made larger for many reasons, one of which was to im. 


٣۶٢ operations efficiency. To realize the full benefit of the larger 


7 


ize, the ship had to be steered as efficiently as possible. Almost 
every Super. Tanker built nowadays is course unstable and the control 
of such a system has been of great concern. We can see easily that 

4 
it is not difficult to design an autopilot if enough rate control is 
Available |5]. 


The basis for the need of Automatic Control as ٣٣٣٣۶۶ 


number of operator controls, minimum rudder orders for course- 
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keeping and automatic course correction without overshoot. But the 
main reason for this type of control is that steering an unstable ship 
by manual control requires a helnmisman to detect.a very s malla: 
of rate of turn and to react to this in time with the proper correction. 

Even though the study of man as a pilot of giant tankers is rel. 
atively young, yet it is evident that any pilot, not withstanding his 
experience and seamanship, has some physiological and psychological 
limitations relevant for the execution of his task. 

Several studies carried out [6 and 7] provided us with some in- 
sight. First a purely physiological problem is a limited capacity for 
Enesperception of motion. Very slow motions cannot be perceived by 
human beings as is the case with the hands ofa clock. To be more 
specific, the smallest yaw velocity human beings can perceive is 
amt I minute of arc per second, The smallest just noticeable 
meecreration or deceleration of the yawmmotion occurs when thie yaw 
velocity is doubled or halved within 5 seconds. Analysis of ship man- 
euverability 0 on full-scale models revealed that tankers 
ےت‎ 100, 000 Ton. move so slowly during maneuvers that many of 
the accompanying velocities and accelerations are not perceptible 
for the man on the bridge. This means that some information which 
is indispensable for the correct execution of the maneuvers is not 


available for the navigator due to his own physiological limitations. 
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The inertia of large tankers also raises a problem of a more 
psychological nature. It is well known to pilots and Captains that 
there is a time lag between a rudder command and the reaction of 
٢1 ۰۳۱1۱ This lag can be easily of tne order of 15 seconds, ra = 
chological research on tracking behaviour revealed that steering 
tasks raise difficulties when the time lag is longer than 4 seconds. 
Also related to inertia is the ability to anticipate motions over long 
time intervals. Even very simple maneuvers take a long time and 
should therefore be initiated a long time beforehand. It is however, 
very difficult to anticipate maneuvers over long time intervals as 
human beings have only a limited capability for extrapolation of slow 
motions, i.e., the navigator can hardly predict what his position 
will be in 10 minutes and therefore is not able to initiate the neces- 
sary correction at the required moment. 

Still another factor is the ability to discriminate among the 
different components of the ship's motion, those due to wind, current 
and steering movements initiated by the pilot. 

nnomsezperiences 1 workmeeswiiitiherconeepr ol Ineiıme ea, 
stent of a ship, conclusions resulted in that the difficulty in steering 
depends largely on this parameter, since for the same statistical 
characteristics, the unstable ship with the small value of absolute 


time constant, changes its course too quickly to be followed by the 


human operator. 


TZ 








II. COURSE-KEEPING 


ar DEFINITION 

Course-Keeping quality is defined as the ability of a ship to keep 
course easily ona seaway. This is not identical to stability on course, 
because the former depends upon the behaviour of the steering devices 
as wellas of the ship, while the latter is a thoroughly passive char- 
acter of the ship. 50 the problem of course-keeping should be treated 
pusmnalyzing-a closed-loopteedback system composed ofa ship, an 
automatic or manual steering contro! unit and steering gear. It is 
possible, however, to show simply that good stability on course 
fetes a pood course-keeping quality, in general, in this torm: 

1。 A ship that is more stable on course has less occasion to 
deviate significantly from its course because the stimulated motion 
decays more quickly even with the rudder amidship. 

2. A ship that is more stable on course responds more quickly 
to steering, so that any BER deviatión 1 9 ۰ ۳ ۱ more 
easily. 

The stability on course is related to the decay of the yaw and 
sway of the ship with a rudder amidship after being disturbed fora 
short while. In the case of an unstable ship this results in a negative 


damping of yaw and sway. Stability on course is also related to quick 
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response to steering. This is of particular importance in considering 


the course-keeping quality of a ship. 


0۶۷۵م 

The dynamic stability on course of an unsteered ship is measured 
I| she number of ship lengths traveled by a stable ship, in the time 
required to reduce an accidental deviation from undisturbed motion 
ina straight line to l/e of its initial value, the rudder remaining in 
the amidship position, this quantity is defined as 1/ | ار‎ : 

A negative value of pj means that the ship is dynamically stable, 
i.e., when a dynamically stable ship moving on straight course is 
disturbed slightly, it settles down ona new straight course that is 
Brose to the original one (Figure 1). The greater the negative mapni- 


Bede of p,, the more rapidly the ship settles on its new course and 


1? 
the closer the new course is to the original one. A positive value of 
p, means that the ship is dynamically unstable, i.e., when a dynam- 
ically unstable ship moving on straight course is disturbed slightly, 
it will go into a steady circuling motion, even though the rudder is 
held amidship. 

The primary purpose of an automatic steering device, just as of 
a human helmsman, is to compensate for disturbances, making the 
ship maintain a prescribed heading and thus making it directionally 


Stable. 
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1 سا۸ رچ و لم‎ DYNAMICALLY 
course  STRABLE > ريم‎ 2 
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“YVSTHELE SHIP 


JJ ( 


eee UNSTABLE HULL 
In this study of Automatic Course-Kecping we will be working 
with a 200, 000 DWT super-tanker of the following characteristics: 
Length 310 meters 
Breadth 47.16 meters 
Draft 18.90 meters 


Steering Quality indices: 


ur = -269. 3 seconds 

seconds‏ 9.3 - را 

T4 = 20.0 seconds 

K 50050434 2d/sec 

Maximum Rudder deflection 30° 

Maximuin Rudder rate 2.32 degrees/second 





D,. AUTOMATIC CONTROL AND COURSE ٔ 1+٦ 

One of the functions of ship control is to maintain a ship's 
heading. In performing this function, a helmsman deflects the rudder 
in a way which will reduce the error between the actual and desired 


heading (Figure 2) 


۷ وج ,جع‎ 
vU. 
e a - 
HEADING 





۷ ACTUAL 
HEADING 


77 AL 


A good helmsman will not only deflect the rudder in response to 
the heading error, but he is also sensitive to the angular velocity of 
the ship and he will ease. off the rudder and apply a little opposite 
E in order to prevent overshooting the desired heading. It 
follows that an automatic control should also be responsive to control 
signals measuring both error and angular velocity. Thus a rudder 
under automatic control might be deflected in accordance with the 


following linear expression: 
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| E Tu 


where 
Se h F helmsman's rudder angle 
e angular velocity of ship 
024 Heading Clio 


The basic action of the aulomatie control or “autopilot” is called 
proportional control which means to give a helm angle proportional 
to the amount of course deviation and scmetimes also proportional 
to the time rate of the deviation. 

In general the complete steering system can be represented as 
in figure 3 


EXTEREAL 
DISTURBANCES 
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CS 


The response of a ship to steering or, the maneuverability of a 
ship is usually described by a set of equations of motion of side- 


drifting velocity and turning-angular rotation: 


I 





4 V u Ay, Aye V 0 E f 
dé ia Qt 422 O b; 


It is however, more convenient to use a single equation of motion 
0515776 immediately a relation between turning angular motion ana 
steering, because control signals in autopiloting relates only to turn- 


ing angular motion. 
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MLS TRANSIER ON DION TING I MN GT 


The general principle of 11563 12۸۱۳۴ hole alone noto ol a 
ship in steering may be summarized as follows: 

1。 The coupling of surging to sway and yaw is ignored because 
of its second-order effect on the latter motions [19]. 

2.  WMydrodynamic forces acting upon the ship are expressed as 
We functions of angle of drift, 1.¢., sway velocity, curvature oF 
mye ahip's path, i-e., yaw angular velocity, and ne of helm. 

- ۱ 

EE ros neral tobees are expressca as linca: ic Ho oO 
the acceleration and angular acceleration of the ship. 

Linear equations of motion of this kind have been introduced in 
a number ofarticles [8, 9, 10 and 17], and have proved promising in 
interpreting and predicting behaviours of ships in steering. In the 


Meration uscd here, they arc written as: 


/ / ( 2 ۲ ٠ e 了 
E on m) A YB = (e meme YE = 
111 ا‎ 
uf) (Ti) 2 7 p a 
— zx EPA 2 (Er) Aa € — Na MB Ad 
Using these equations we can obtain the yawing (© (t)) and the 
swaying ($ (t)) in response to any movement of ihe rudder ( á (00 
Now limiting our interest to the steering-to-yaw response which, 


in general, is naturally of major importance, we introduce the 
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transfer function in steering. Taking the LaPlace transform of both 


sides oí te equations Ll و‎ ١٢ ٢/٢ 


-(( Vir) (n moy) PU) “lt Sr) (wie ティ ys] Bos) 
(Y y) (mame =) OCS) = Fg dcs) 


(ZN er Te) Sy Al (Tla tTa) S (HIE ] 0) 
BEE 


where f (0) and 6 (0) denote f and ei respectively at t=0, 
namely at ue beginning of a rudder movement, and where 
P (3) is the LLaPlace transform of f? 
= (S) 1s he Caria C 0: ②③ 
E (e) وز‎ ۱١ را‎ transiorim of á 
laking into account the Burn relation at t=0 
Su) (Liz e Jas) Ec) + (Yo) 04 8 to) - Mg T 
` の 
we can eliminate f(s) and P (O) from the transformed equations [19], 
mu Teld a single equation that describes the steerinp-to-yaw re- 


1 بجی‎ 
Sponse, "that is: 


geen Kgl, ل‎ )5( E 人 AGT] Oto) +T, Ta Gre) 
(Uis) (o T5) (^* 05s Cl tT s) 


MISZ 


jo . ۱ ۱ ۱ : 

Since later in this study we are going to be dealing with auto- 
matic piloting, and the control signal under this condition relates 
only to turning angular motion. 


60 





where 





cl yf ww 
K =(%) Hr Ys + Ya Ys 


/ / / 
Ya Wr 一 (mi m= Yy) Ai م‎ 
er! / / ーー/ BS / 
۳ Tom A ( FA e My) rv + (Liz ter/ قر‎ 
و‎ 7 7 
Va Ne ー( 291 Pts -Yr)^ m 


/ 
Cm hy] Cue Jii) 


T, Ta = (2/1) 


/ 2" 
Ys My (m + KA, m ) 4 


/ 
( P Lat パケ 


Br = Cz ) DE po AN A/ E 
(3 Ge E /rz کے‎ 
The first term of the right side of equation III-2 corresponds to a 
motion excited by steering and the second to another resulting from 
the initial motion at t-0, the latter disappears if the ship is running 


straight at t=0 without any yaw or sway. Thus the rational function: 
€ cs) K (1 Tas) 
ENS ) ۱+۲, 5( )۱+۲5( 





describes the response behaviour of a sbip in terms of the LaPlace 
Beamstorm. We may call it the transfer function of the ship in steering. 
Retransforming equation III-2, we obtain a single differential 
equation that describes ship motion just as does the original equation 
of motion (Equations III- 1) as far as the مد‎ response is 


concerned, that is: 


d, fa 〆 の 2 Ern) 483 e MIAK d 
LB = LUNG At 





ER 





better known as Nomoto's Equation, this describes directly the yaw 
mesponse of a ship. The four coetiicvente dc, ie ول‎ and T3, which 
are composed of the coefficients of the original equation of motion, 
constitute a set of characteristic figures representing the response 
behaviour of the ship, they are called the steering quality indices. 
۲۰۰۰ ٠ LABILITY CRITERION OF AUTOPILOTING BY ROCHENBURE ERS 

METHOD 

A control loop describing autopiloting of a ship is illustrated in 
figure 3 in the form of a block diagram. Since the stability criterion 
of Kochenburger [11] is based upon whether sinusoidal signals grow or 
decay in circulating through the loop, it is necessary first to obtain 
the response of each element composing the system, to sinusoidal 
signals. 

The response of a ship to a sinusoidal signal (that is, in this case, 
to put a rudder sinusoidally to both sides with a certain frequency(w) 


Day be determined through the equation of Nomoto: 
の IPM 1 Son wt 
OC AOI [we + の 227] 


where 


Aug | KA rin) 
Cl*cwT, ) ) ١٥» زبس‎ 


の [し 大 Aro K EN RN CU 
CLAN) AUS 





Ze 





and where 7 is the amplitude of sinusoidal steering. A(W) is called 
an amplitude ratio and の (w) a phase difference. Both of them are 
functions of frequency w only, as a common feature of lincar systems. 
Next, the response of an electrohydraulic steering gear, which is 
widely used for most present ships, may be described by the following 
equation: f 
n グ -f کے‎ = d 1 
uere lis a time constant of the steering gear, 
Then we get a description of the response of a steering gear toa 


sinusoidal JÛ’, as follows: 
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BR. Des "UN CO CE 


A ー Ai, Cw) Dec [wt + の (0) | 


where 


A. (0) = / 


۱ م7 دب + / 


の (w) = Ara ou A 


Finally, we get easily the response character of a compass as 
follows, because it may be considered a simple integrating element 


transforming © into © : 
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2 ep 


OD = Oo Hn vt 
ند سس ۸ - م8‎ A 2. 


ا 


where 
7 
d Ww) 一 ل‎ = 
/ و‎ 77 
Gau: Ze Dodo eS 


Since the transmitting character through all linear elements has 
been thus obtained, if à similar character of the remaining element, 
viz, an autopilot, is defined, we can judge whether sinusoidal signals 
SIO or decay in cireulating through the control ioo It is impossible, 
however. 2 describe the response of an autopilot by any linear differen- 
tial equation and then to obtain its response character in the foregoing 
manner, if considering such a discontinuous element as a weather- 
adjust mechanism, which is used to modified proportional control in 
order to avoid frequent steering in practical application for rough 
seas. 

Fortunately, however, the response itself of an autopilot with a 
weather-adjust mechanism to a sinusoidal input signal (that is, in this 
case a course deviation) may be obtained easily by expanding the re- 
sponse into a Fourier series with a fundamental frequency that is the 


same of the input signal, as follows: 


T ۰ 
Ó = C, K (U, das の C FA, co wt * 3 pev 2 نه‎ «4 Coo 2e...) 
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where 0 is the proportionality constant connecting a 
1 course deviation to a helm angle to be called 
OF, 
2 is the amplitude of course deviation, and 
a,a,a, anda are the coefficients of Fourier series which 
4 

may be obtained through the usual procedure 
of Fourier expansion if the form of A” is given. 

sidering here that all linear elements are much more 6 111- 
tive to a signal with the higher frequency, we can neglect all the higher 
frequency terms. ‘This is the approximation of Kochenburger, and its 
validity depends upon how much the higher frequency signal decays 
through the linear elements and also how much higher frequency come 
ponents are included in the original の ; viz, how much the original S 
resembles a pure sinusoidal form. Inthe present case, Kochenburger's 
approximation may be fairly valid because a ship is quite insensitive to 
a high-frequency steering because of her large inertia, and also the 
basic action of an autopilot, even with a weather-adjust mechanism, is 


proportional control that produces a sinusoidal output in response to 


a sinusoidal input. 


B. FIRST ATTEMPT TO STABILIZE THE SYSTEM 

The autopilot, the steering engine and the ship with its rudder all 
form different components of a closed loop Sale each component 
characterized by its transfer function or the complex ratio of output to 
input. The theory for such control systems and their stability has been 


developed in electric network engineering, and it is natural that the 
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dynamic problem lends itself to studies 7 217 191 9 ١٢ Dee 
each component is represented by its equivalent electric circuit. The 
ability of the closed loop system may be judged from the total open 
loop response recorded at several frequencies, without a knowledge 
of the individual transfer functions. 

The conditions of directional stability in automatic steering along 
a fixed course had been discussed by Minorsky [12] by means of the 
technique of added derivatives, applied to a simplified one-degree 
of freedom oscillation and including severaltypes of position and rate 
control. In 1946 Davidson and Schiff [13] took a large step towards a 
better o رکب‎ of the interrelation between the performance of 
a ship on a straight course and in turning, pointing out the nonlineari- 
ties in the behaviour of the unstable ship and including an interesting 
treatment of the transients when entering a steady turn. They also 
used linear theory for establishing a formula for the radius of the 
Maning circle. 

The work of 11 and Schiff did much to stimulate other 
Enos. Among many others, a paper by Williams [14] on initial 
stage motion and a report by Schiff and Gimprich [15], who studied an 
automatic control system where the rudder angle called for is pro- 
portional to a combination of heading deviation and rate of change of 
heading, and which has a behaviour with a close resemblance to the 


automatic pilots used in practice. 
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Most modern steering engines are designed to move the rudder 
with an essentially constant speed, the rudder turning at that speed 
as long as a control signal is transmitted to the stecring engine. A 
gyro pilot may be used to switch on this signal, calling for a correct- 
ing rudder at a certain small deviation from the desired heading, 
whereas a contact on the rudder may stop the motion at a suitable 
angle, the rudder remains in this position until the ship has swung 
Over to the other side of its course and the rudder is then reversed. 
More often a follow up mechanism assures a "proportional control"! 
ei the Eder: Due to the finite rudder speed, these systems may be 
Siren E citing, however, and excessive oscillations may be built up. 

In order to overcome these difficulties the automatic pilot must 
be made to anticipate the motion of the ship, much in the way an 
experienced helmsman gives an auxiliary rudder. In practice this is 
accomplished by means of a feedback of rudder motion to the heading 
mie detector, as in many commercial applications, sometimes also 
by adding some kind of rate of change of heading control, similar to 
the Mich papclcobupeneppodlsubrdjarinoc'depPbh control sySteims..) in 
effect, both these methods correspond to a "proportional plus rate 
control;'' the first one often incorporates a, non-linear character to 
the system [16]. 

With these concepts in mind we are going to try to stabilize our 
unstable system, first by assuming that the autopilot is just a gain G, 


we represent the block diagram in figure 4 
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FIGURE 4 


The open loop transfer function for this svstem is: 
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Since there is one pole in the right hand plane, given by TE there 
is no sense in ۶ 13۳ علا‎ 250106 DIOL {GS 171-۸116 18 ٦ 

Il.  Root-Locus Analysis 


Gettinp the characteristic equation and arranging it in the 


proper form we obtain: 
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Since we have four poles there are four separate loci. 


There are three asymptotes due to 


n#of poles = 4 
n-q = 3 
مم و‎ of zeros = 1 | 
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From a sketch made for this analysis we get the cross point 
at the jw axis = .025. So Gat the crossing point: 
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To find a better set of mambiers we reter to the Routh 


criterion: : 
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replacing the values of G, and after al] the mathematics we obtain: 
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values that agree with those approximated trem the sketch. 

In the figure 5, computer output of the subroutine Root Locus, 
we obtain the open loop characteristic of this systern. As we can see 
the system can be considered as marginally stable due primarily to the 
long transient that can be predicted from the root locus graph. 

Using the DSL package program we simulate the system, 
first the system represented by the block diagram of figure 4, and 
after that the same system with the concept of proportional plus 


rate control, ' using rate of turn as feedback, as represented in 
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In both cases, the disturbance was simulated as a step of 
amplitude 1, values of G were selected between 12 and 36, which are 
the ones with better time constant, copy of both computer programs 
is given at the end of the thesis. 

Figures 7-12 are computer outputs of the first condition, 
meen, System of figure 4, we can see as expected, the Jong time 
required to settle down due to the small time constant Of the G values, 
but even so we can see in figure 9 for the value of G of 24.2, com- 
paratively a short time for reaching steady state. 

Figures 13-16 are computer Outputs of the second conditions 
رہل‎ Eon of figure 6, we get better damped results, with a sig- 
nificant decrement in the time to reach steady state as expected. 
Also in figure 15 we can see that comparatively this value of G pives 
us the best results. So we can state that the value of G = 24.2 is a 


good value to stabilize the system and perbaps the best under this 


conditions. 
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IV, Besic COURSE-KEEPINGREIUEE MEN 


It is well known that the course-keeping quality and the turning 
quality of a ship are contradictory to each other, and that it is im- 
possible to improve both of them without increasing rudder area [4]. 

The success of a rudder, designed for à particular ship, is 
measured by the degree to which the ship achieves a desired or anti- 
cipated course-keeping and course-changing ability. Since most ships 
spend most of their operating time moving ahead, it is as i that 
Maneuvering performance is specified only for the ahead direction. 

Course-Keeping ability may be specified quantitatively in terms 
Beatie stability index, of the characteristics of the Dieudonne Spiral 
Maneuver (see Appendix B), or of the range of rudder angles used 
to maintain a straight course. e In reference [18], itis suggested 
that an attempt be made to design all ships for a stability index of 
ECHO or less, but itus recognized that this may not be practicable for 
all ships. 

A major reason ior the specification of a level of control-íixed 
stability for all ships is the fact that excessive controls.fixed 


instability leads to excessive use of the rudder to maintain a straight 


CET‏ : ت2 
This measure of Course-Keeping ability is frequently used in‏ 
maneuvering in restricted waterways or in rough seas.‏ 


46 





course. This in turn leads to increased wear on the rudder system, 
decreased ship speed, or increased fuel consumption. Controls- 
fixed instability also leads to increased difficulty in navigation in 
restricted waterways and in following seas and in avoiding collision 


with passing ships. 


ESS UNDER STEADY CONDITIONS 
1, Calm Water 

Among the factors important for the turning motion of ships, 
there are transverse forces and moments acting upon ship's body 
itself. These may be approximately looked upon as motion in infinite 
fluid where image of ship is considered in case when V is not so large 
and the effect of waves by the ship itself with Froude Number below 
about 0.2 or 0.7 ۷ L/2d [16] is negligible or when motion of ship 
is slow, and are the functions of V, f E の . This will be sufficienti 
for turning motion of ordinary ship. 

The disturbances acting upon the system can be classified 
according to their nene on the behaviour of the system as follows: 

a. Pis that cause deviations from the set course. 

b. Disturbances which affect the steering characteristics of 
the ship. 

Wind and waves belong to the first group, to the second class 
of disturbances belong the loading of the ship, the depth of the water, 


CIC. 
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In considering the calm water condition we arc assumüng no 
disturbance at all, so, under these conditions we expect the ship to 
maintain its heading, without any action of the autopilot. 

2, Steady Wind 

۱۳۱ ۰۵ the wind pressure, the chip لا‎ ۵ 2 6۱ upon bya 
lift and drag which will cause the ship to drift at a drifting angle. By 
the drifting, the hull together with the rudder will create a certain 
hydrodynamic lift, drag and moment until the windward component of 
the hydrodynamic force balances with the wind drag. At this condition, 
the total summation of forces including ship's thrust is zero. How- 
ever, the summation of the moment including moment due to the rudder 
at maximum helm angle will not necessarily be zero. 

If the summation of the moment is not zero, and is ina direc- 
Bon tO increase, the ship will have a tendency to turn windward. If 
the moment is in an inverse direction the ship will turn leeward. In 
both cases, the ship is supposed to be uncontrollable, If the summa- 
tion of the moment = zero within the range where ruddcr angle does 
not E eed 30 degrees, the ship can keep its heading so that it is 


supposed to be controllable. 


B. ACCURACY NEED FOR COURSE-KEEPING 
۳ Loop Gain 
The most important parameter to be considered in considering 
automatic control is the Loop Gain, and its value depends directly on 


the transient response of the system. 
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For linear systems, the most commonly used correlation 
between frequency response and transient response is that correlation 
which exists between the height of the resonance peak Mog and the 
height of the peak overshoot of the step response Ss for a second- 
order system. A curve relating these is easily calculated [20]. 
mcn, for a known second-order frequency response, the curve is 
entered with Mpw and the transient peak overshoot Mpt is predicted 
Sly. When the system is known to be third-, fourth-, or higher- 
Order, there is no such readily available correlation, but the second- 
order correlation may be used as an estimate. This gives da very 
accurate estimate if the higher-order response is dominated by one 
pair of complex roots, a less accurate estimate if dominance is not 
assured. 

We know that in order to make a system stable within some 
degree of accuracy we need to have complex roots, and a certain 
value of gain. In our study where we are operating with a giant of 
200,000 Tons., we 5 aora 0 1376 ٥ 1ک‎ 1126 Orna للا‎ 
Ge , 1 tO 2ugegreesset headamp error, withoutteensequences, So we 
ameo, 11 required, to have real roots governing our system. 

2. Loop Type Number 

Control systems, in general, are required to have certain 
operating characteristics which are determined quantitatively by the 
aP eiie control problem, but which may be listed qualitatively the 


ones we are interested in as follows: 


A) 





a. The system must be accurate in steady stat, 

b. Ihe 527610617171011 be stable, 

c. The system must regulate against disturbances. 

Each problem onsfeedback eonlrol has a different set or vers 
formance requirements that must be satisfied, but the requirement 
HEUS common to virtually all problems is a need for accuracy. 
What is meant by accuracy depends on the specific physical applica- 
tion, and no convenient broad definition is available, but the scope of 
accuracy requirements is readily illustrated. The class of control 
system we are interested in is the class of positioning system or 
servomechanisms. These systems are designed to change the output 
quantity as commanded by an input signal, and in addition are required 
to act as regulators in the presence of output disturbances. The 
primary consideration is usually that of keeping the error (difference 
between commanded output and actual output) less than a specified 
amount when the system is in steady state with no load disturbances. 
If static load disturbances are anticipated a separate specification 
for NE us error is usually given, and for a suddenly applied disturb- 
ance a maximum instantaneous error may be specified. 

For positioning systems subject to command input signals, 
there are also specifications as to the permissible are of the 
transient response, but these specifications usually refer to the 


permissible time duration of the transient and to the permissible 
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nature of any oscillations during the transient period. They seldom 
refer directly to the accuracy of the system, but the constraints 
which they place on the system design may make it very difficult to 
satisfy the accuracy specifications. 

Systems classified as positioning systems may be subjected 
to various types of input commands. These may be well defined and 
representable as simple mathematical functions or they may be rather 
eemplex functions of time. for well defined inputs the system usually 
(not always) reaches a finite steady state condition so that the accu- 
racy of the system is readily evaluated. For more complex inputs 
Bes stem may not reach a clearly defined steady state condition, 
so that the accuracy must be evaluated in terms of bounds and a 


numerical value for "steady state error" 


may not ۰ 
In this study emphasis is placed on system analysis and 


design for cases where the accuracy problem is defined in terms of 


deterministic input signals. 
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Having been able to obtain the best value for the gain of the 
autopilot to stabilize the unstable model, we see that the transient is 
still too large and has overshoots (due primarily to the hydrodynamic 
coefficients of the ship, which are given, as we may recall by the 
steering quality indices). 

With the concepts of chapter IV we are going to approach the 
design of the compensation using the Root- Locus method instead of 
the Bode Plot, since having a pole in the right hand plane the latter 


is hard to use. 


۱ نز ی‎ GYRO COMPENSATOR 

Because of the shape of the Root Locus diagram we need to 
improve the transient response. By using a Rate Gyro we introduce 
an additional zero at 1/k, and this relocates the roots on the Root 
Locus plot ata deser Balin level. Kate feedback 1s a Very CONO 
means of increasing equivalent fiscous damping and of thus improving 
system transient response. Figures 19 and 20 are a block diagram 
representation oí heading feedback, the original circuit and the 
equivalent respectively. 

The new characteristic equation is 
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we obtain the 4th order polynomials to use in the Root-Locus 


analysis: 
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Figures 21-26 are computer outputs of the subroutine Root-Locus. 
An sive the different system behaviour for the values of k selected. 
As we can see from figures 24-26 in order to have a good time constant 
we need a large value of G, and from figures 21.23 we can select 


very small values of G, which will give us real roots. 


fee AD TYPE FILTER COMPENSATOR 


Another type of compensation uses a filter with transfer function 
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of the form: = "د‎ 
sg を 7” 
AA سا يي‎ 
© + م‎ 


having high pass characteristics. We are going to relocate, also, 


Bie TOOLS Ol the Noot-Uocus diagram.  Thisfiltemus just a 8ی‎ 
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transfer function in the direct path between input and output, on our 


system. We are going to introduce one additional pole and zero 


Figure 27 shows the block diagram with the filter in the system 
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Figures 28-33 are computer outputs of the subroutine oo Locus 
iimey pivesthe different system 06٢3۷11 ۱۱۴١ ٤ ۶ ٦ 
selected. We can see from figures 31-33 that in order to have a 
good time constant we need a large value of G, while from figures 
28-30 we can select very small values of G, which will give us real 
roots. 

Having obtained the computer outputs for both types of compensa- 
tion, we can observe à similarity in the system behaviour for same 
root location. So we can predict the same transient for step inputs. 
However, we need to find an optimum value of k in the case of the 
ire Gyro and of z for the Filter, in order to expand the little loop 
in figures 21-23 and 28-30. This optimum is the best value possible 
to obtain and it is a compromise between a high gain and a good rudder 
operation, bearing in mind that for this to occur we must consider 
maximum deflection and rate for the rudder of this ship. Since fora 
G@uick response we need a high gain, but this could demand an exces- 
sive operation of the rudder and this could be bad for the following 
E... 

1l. ‘To get quick response we need faster rudder angle and a 
۶۰. )6 rudder angle. 

2. Faster responses may cause accelerations that are too 


sharp for personnel on ship. 
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So the next step will be to try to find the optimum value mene 
tioned above. Right off hand we can consider {he system benaviDums 
erioures 25, 26, 32 and 33, are very much similar to our original 
System without compensation, so we can disregard those four systems, 
and now we have four different systems for the Rate Gyro and Filter 
Esrmopensators respectively. 

We Se going to simulate a disturbance that produces a rate of 
turn of . 2 degrees/second, a value that has been found by experience 
can easily be reproduced on this tanker, requiring a correcting rudder 
angle of 10 to 20 degrees. 

Maith this external disturbance represented as a Step, we examined 
our eight compensated systems. By trial and error method we check کل سا‎ 
for the optimum value required and we end up witha I/k =. 041 for. 
Mee cte Gyro, anda z = . 041 for the Filter, which in essence are 
Re same values we have for the case of 1/k = .04 and z = .04 for the 
¡Mate Gyro and Filter respectively, so now remain these two systems 
to be studied for the bed gain. Inall the eight initial cases we ob- 
served a common characteristic, with minor variations, that fora 
small variation of heading, a sudden change in rudder occurs which 
causes a very quick response and we may suspect an excess in both 
۱۳۱۳ 6۶ parameters. 

Figures 34-37 are computer outputs for the system compensated 


with the filter, we have chosen values of gain from the root-locus 
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graph, for four different damping coefficients, 9.42, 15.87, 23.67 
and 3.5 respectively. We can observe from the first three figures 
that a quick response is obtained and steady state is obtained ina 
ur short time, but this could exceed theo restrictions of the ۲ 

In figure 37 we see a moderate response and we do not reach the 
steady state at the end of the simulated time, which is only 5 minutes, 
as expected for a ship of characteristics such as this. We need to 
verify if there is any violation of the rudder restrictions. 

Matin the help of computer simulation, a program im Dol, which 
is given at the end of the thesis, is used to simulate limiters for both 
parameters. In order to introduce the limiting value of Rudder Rate, 
we need to replace, in the block diagram, the steering Gear transfer 


lunet!on fora circuit equivalent, given in figure 38: 





FIGURES 
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DIXI Um 
3 B 


With the program above mentioned and the values of gain selected 


before we obtain the figures 39-42, computer outputs, where we see 
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in the first figures that even if the system is stabilized, there is a 
portion of the response where we exceed the rudder deflection, while 
in figure 42 we see a smooth settle down without any clipping. Will 
be a good idea to try values of gain within this small range of gain 
values. 

Figures 43-48 are computer outputs for the system above men- 
۳/۵۰ with values of gain of 6., 5.3, 4.6, 3.5, 2.5 and 25 Wecan 
observe in figure 44 for a gain value of 5.3 we obtain a good settling 
time with the smaller steady state error of only .8 of degree out of 
course which is excellent, but we have some rudder action even 
after settling time. In figure 45 for a gain value of 4.6 we obtaina 
Sacady state error of 1.03 degree out of course, with a settling time 
of about 2 more minutes than for the former case, but we get smaller 
rudder deflections. So we can conclude that this value of 4. 6 will be 

E ا‎ value for this system. 

S T cing the Rate Gyro having a zero at 1/k = .04, with 
the same disturbance, we obtain the computer outputs of figures 49-52, 
دنا‎ 2165 of gain selected from the Root-Locus graph, for four dif- 
ferent damping coefficient, 22.26, 67.28, 35.29 and 3.5 respectively, 
Epobservesdlirsospthessamerresultscas imme cae ec Ol Mc HICET, l.C, 
the first three figures evidently violate the rudder restrictions and 
only in figure 51 for a gain value of 3.5 gives us a good system 


operating under real conditions. Figures 53-58 are computer outputs 
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for this system, for gain values 6i 0920500 05 4 45 8 ٢ اکا‎ 
following the same considerations as with the Filter. We can see 
from these figures that in all des the system has increased its 
transient, being the best value, 6.2 as observed in figure 53. In the 
event of increasing the gain to 6.7 we start obtaining clipping in the 
rudder action, so this value of gain is chosen as the best suited for 
this type of system. 

From comparison between these two systems we can select the 
one with Filter compensation, first, because of the best final results, 
i.e., less amount of rudder required, better settling time, and 
secondly E the reliability of the system, since the Rate Gyro could 


nube sultable for marine use due to this overhaul interval. 


\ 
تست د هه ^ ^ | SS‏ 
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VI. CONCLUSIONS | 


۰۰۰٠٠۰۱ ٦۹۷٢(۵ 

An optimal value of gain can always be obtained using the Root- 
Hocus method for any given plant. Both negative and positive feed- 
back can be used to achieve the desired results but only negative 
position feedback should be used to obtain a stable system. 

With the concepts of Chapter IV, the requirements for a basic 
Course-Keeping, and our best value for the system selected as the 
best, we verify with a low order disturbance represented by an initial 
condition of © , that the actual heading in the steady state is equal 
to the reference heading or which is the same the steady state error 
is equal to zero. We already have observed that for a load disturbance 
represented by the step input we obtain a small steady state error. 

So we have reasons to believe that we are working with a type zero 
system and, ifitis so, do we need to change the loop type number? 

By doing an an ofthe feedback loop we can easily verify 
that indeed we have a type zero system, and using the final value 
theorem for this linear system, the evaluation of steady state error 
for the deterministic input is easily a 0 In this case 
where steady state error is finite its magnitude is determined by the 


GEeCciprocal of the gain transfer fumetion for therdinee: path between 
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input and output, this magnitude is called an error constant, also 
called the DC gain or zero frequency gain, because it is the numerical 
value obtained for the transfer function by deleting the s factors and 
substituting s - 0. Following this reasoning and with the final results 
obtained we see that this meets our requirements for the stable 
operation of the system with the degree of accuracy that can be con- 
sidered as satisfactory. 

The purpose of this thesis, which was to study the basis of the 
Course-Keeping with Automatic Control, is reached at this stage. 
It is good to point out that further studies could be made by studying 


the sea effects, how automatic steering and the yawing of ships is 


affected in rough seas. 
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APPEND A 


LIST OF S TMBOLTS 


draught of a ship, mean of fore and aft if no remark. 


moment of inertia of a ship about a vertical axis through 
Mer centre of gravity. 


ZE I 


moment Cr inertia Of additional mass about a vertical 
through the centre of lateral additional mass. 


T ا ان‎ 
zx کم‎ 2 
[ 1180 Dip, between perpendiculars if no remark. 
mass of a ship. 
E 6 

a 
longitudinal additional mass. 
p ۶َ ٤ 

£ 1 
Jateral additional mass. 
LE 

e 


hydrodynamic moment about a vertical axis through the C.G. 
Ola ship, positive for standard turning moment, 


۸۷/ + t 


du A 
IND y 


0 2 が 


yaw angular velocity, positive to starboard, identical with © . 


ae 





په 


< 
c 


D ہہ‎ D 


1۸۲ ۷ ۰۰۶۶ یی‎ Ende 
turning radius, steady or instantaneous 
ship speed in meters/second 


lateral component of hydrodynamic force, positive to star- 
board. 


t 
تس و‎ 
í 
OI don 
اس دال‎ 
/ 
۵2 däs 
angle of drift, positive to port. 
angle of helm, positive to starboard. 


turning angular velocity, positive to starboard, identical to r. 
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THE 0171000000۶5 ۶+٣ 


The spiral maneuver consists of the following: 

1. The ship is "steadied" on a straight course ata preselected 
speed and held on this course and speed for about 1 minute. Oncea 
steady speed is established, the power plant controls are not manip- 
ulated for the duration of the maneuver. 

2, After about 1 minute, the rudder is turned to an angle, A : 
Si 5150111 15 degrees and held until the rate of change of yaw angle 
maintains a constant value for about l minute. 

5. the rudder angle is then decreased by a small amount 
(about 5 degrees) and held fixed again until a new value of e is 
achieved and is constant for several minutes. 

4. The foregoing procedure is repeated for different rudder 
angles changed by small increments from, say, large starboard 
values to large port values and back again to large starboard values. 

The numerical measures obtained from the preceding spiral 
maneuver are the steady yawing rates as a function of rudder angle. 
A plot of these values is indicative of the stability characteristics 
wea slip. Por example, if the plot is 3 single line going from star- 
board rudder to port and back again, as shown for ship A in figure l1- B, 


the ship possesses controls-fixed, straight-line stability; that is, it 
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1002 


has a negative stability index. If, however, the plot consists of two 
branches joined together to torma "hysteresis" loop, as shown in 
ship B of figure 1-B, the ship is unstable; that is, it has a positive 
stability index. In addition, the height and width of the loop are 
numerical measures of the degree of instability; the larger the loop 
the more unstable the ship. The slope of the yaw-rate curve at zero 


rudder angle is a measure of the degree of stability or instability. 
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